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Role  of  lee  waves  in  the  formation  of  solid  polar  stratospheric 
clouds:  Case  studies  from  February  1997 
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Abstract.  Recent  theories  of  solid  polar  stratospheric  clouds  (PSCs)  formation  have  shown  that 
particles  could  remain  liquid  down  to  3  K  or  4  K  below  the  ice  frost  point.  Such  temperatures  are 
rarely  reached  in  the  Arctic  stratosphere  at  synoptic  scale,  but  nevertheless,  solid  PSCs  are 
frequently  observed.  Mesoscale  processes  such  as  mountain-induced  gravity  waves  could  be 
responsible  for  their  formation.  In  this  paper,  a  microphysical-chemical  Lagrangian  model 
(MiPLaSMO)  and  a  mountain  wave  model  (NRL/MWFM)  are  used  to  interpret  balloon-borne 
measurements  made  by  an  optical  particle  counter  (OPC)  and  by  the  Absorption  par  Minoritaires 
Ozone  etNO*  (AMON)  instrument  above  Kiruna  on  February  25  and  26,  1997,  respectively.  The 
model  results  show  good  agreement  with  the  particle  size  distributions  obtained  by  the  OPC  in  a 
layer  of  large  particles,  and  allow  us  to  interpret  this  layer  as  an  evaporating  mesoscale  type  la 
PSC  (nitric  acid  trihydrate)  mixed  with  liquid  particles.  The  detection  of  a  layer  of  solid  particles 
by  AMON  is  also  qualitatively  reproduced  by  the  model  and  is  interpreted  to  be  frozen  sulfate  acid 
aerosols  (SAT).  In  this  situation,  the  impact  of  mountain  waves  on  chlorine  activation  is  studied. 
It  appears  that  mesoscale  perturbations  amplify  significantly  the  amount  of  computed  CIO,  as 
compared  to  synoptic  runs.  Moreover,  MiPLaSMO  chemical  results  concerning  HN03  and  HC1 
agree  with  measurements  made  by  the  Limb  Profile  Monitor  of  the  Atmosphere  (LPMA) 
instrument  on  February  26  at  a  very  close  location  to  AMON,  and  explain  part  of  the  differences 
between  LPMA  measurement  and  Reactive  Processes  Ruling  the  Ozone  Budget  in  the 
Stratosphere  (REPROBUS)  model  outputs. 


1.  Introduction 

In  polar  regions,  heterogeneous  reactions  which  occur  on 
stratospheric  aerosols  and  on  polar  stratospheric  clouds 
(PSCs)  are  a  crucial  step  in  chlorine  activation  responsible  for 
ozone  loss  [Solomon  et  al ,  1986].  The  efficiency  of  those 
processes  depends  on  three  parameters  strongly  linked  to  the 
temperature:  particle  type,  surface  area,  and  chemical 
composition  in  the  case  of  liquid  particles.  Thus  the  way  in 
which  PSCs  form  is  of  primary  importance.  Our  knowledge  of 
PSCs  formation  has  evolved  over  the  past  few  years, 
particularly  concerning  solid  particle  formation.  Five  types  of 
particles  are  usually  considered:  liquid  binary  aerosols, 
sulfuric  acid  tetrahydrate  (SAT)  particles,  type  la  PSC  made  of 
nitric  acid  trihydrate  (NAT),  liquid  type  lb  PSC 
(H2SO4/HNO3/H2O)  also  named  supercooled  ternary  solution 
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(STS),  and  type  II  PSC  made  of  ice  particles.  Some  pathways  of 
solid  PSCs  formation  remain  uncertain,  especially  concerning 
PSC  la  [Tolbert,  1996],  Tabazadeh  et  al.  [1994]  and 
MacKenzie  et  al.  [1995]  suggest  that  PSC  la  can  be  formed  by 
freezing  of  STS.  Zhang  et  al.  [1996]  show  that  NAT  could  form 
by  nucleation  on  preactivated  SAT.  Tabazadeh  and  Toon 
[1996]  suggest  that  solid  type  PSC  could  form  through  a 
metastable  dilute  HNO3/H2O  solid  phase.  Other  works  point 
out  that  PSC  la  could  be  composed  first  of  metastable  nitric 
acid  dihydrate  (NAD)  particles  which  could  form  in  a  strong 
warming  event  [7sras  et  al.,  1997],  and  convert  later  into  NAT 
[Middlebrook  et  al ,  1992].  More  recent  studies  suggest  that 
NAT  could  form  by  nucleation  on  ice  nuclei  [Biermann  et  al., 
1998],  Type  II  PSC  would  therefore  be  the  first  step  in  the 
formation  of  solid  PSCs  [Koop  et  al.,  1997].  SAT  particle 
formation  was  studied  from  binary  solution  freezing 
experiments,  although  this  mechanism  does  not  occur  under 
polar  stratospheric  conditions  [Clapp  et  al ,  1997;  Carleton 
et  al.,  1997],  since  the  droplets  are  STS  at  temperatures  at 
which  they  are  supposed  to  freeze.  SAT  formation  mechanism 
is  therefore  still  unclear.  However,  SAT  particles  could  play 
an  important  role  in  PSCs  cycle  formation:  Koop  and  Carslaw 
[1996]  have  shown  that  SAT  deliquescence  could  convert 
SAT  aerosols  into  STS  droplets,  which  could  freeze  later  to 
produce  solid  PSCs  [Iraci  et  al.,  1998]. 

Our  knowledge  about  the  temperature  at  which  solid 
particles  are  formed  has  also  developed  over  the  last  few  years. 
Koop  et  al.  [1995]  have  shown  that  liquid  particles  would 
not  freeze  above  the  ice  frost  point  (hereafter  Tux).  Koop  et  al. 
[1997]  suggest  that  liquid  particles  could  remain  liquid  down 
to  3  K  below  T, ce,  while  Carslaw  et  al.  [1998a]  suggest  4  K 
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below  Tlct.  This  imposes  very  restrictive  conditions  for  solid 
particle  formation,  since  at  synoptic  scale,  such  temperatures 
are  not  usually  reached  in  the  Arctic  stratosphere.  However, 
mesoscale  processes  such  as  mountain-induced  gravity  waves 
can  generate  vertical  air  motions  that  are  large  enough  to 
adiabatically  cool  the  air  parcels  to  3-4  K  below  Tice  [Cars law 
et  al.,  1998a],  The  role  of  mountain-induced  waves  (lee  waves) 
on  PSC  microphysics  and  chlorine  activation  has  been 
investigated  recently  [Meilinger  et  al.,  1995;  Tsias  et  al., 
1997;  Carslaw  et  al.  1998b,  1999], 

In  this  context,  the  aim  of  this  paper  is  to  study  two  cases  of 
particle  observation  in  February  1997  over  Kiruna  (67.9°N/ 
22.1°E),  Sweden,  which  cannot  be  explained  by  synoptic 
scale  temperature  histories.  On  February  25,  1997,  a  balloon- 
borne  optical  particle  counter  (OPC)  [ Deshler  et  al.,  1993] 
measured  a  layer  of  large  particles,  while  the  balloon-borne 
Absorption  par  Minoritaires  Ozone  et  NO*  (AMON)  [Renard 
et  al.,  1996]  instrument  detected  a  layer  of  solid  particles  on 
February  26.  The  Microphysical  and  Photochemical  Lagran- 
gian  Stratospheric  Model  of  Ozone  (MiPLaSMO),  which  has 
been  developed  to  study  the  chemical  composition  of  the 
stratosphere  (taking  into  account  a  detailed  PSC  model, 
stratospheric  chemistry,  and  heterogeneous  reactions  on 
PSCs),  is  used  to  tentatively  explain  the  measurements.  The 
orographic  gravity  wave  model  of  Bacmeister  et  al.  [1994]  is 
also  used  in  this  work  to  investigate  the  role  of  mountain 
waves  in  the  formation  of  these  disturbed  layers. 

The  models  used  here  are  described  in  detail  in  section  2, 
then  the  instruments  OPC,  AMON,  and  Limb  Profile  Monitor 
of  the  Atmosphere  (LPMA)  are  presented  in  section  3.  For  the 
two  cases  considered  in  this  paper,  the  role  of  lee  waves  in  the 
formation  of  disturbed  layers  is  investigated  in  sections  4  and 
5.  In  the  case  of  February  26,  1997,  the  role  of  mountain  waves 
in  chlorine  activation  is  discussed  by  comparing  MiPLaSMO 
chemical  outputs  with  Reactive  Processes  Ruling  the  Ozone 
Budget  in  the  Stratosphere  (REPROBUS)  Chemistry  Trans¬ 
port  Model  (CTM)  calculations  and  LPMA  measurements  of 
HN03  and  HC1  [Pay an  et  al.,  1999]. 

2.  Model  Description 

2.1.  MiPLaSMO 

MiPLaSMO  is  a  numerical  model,  which  was  developed  to 
study  in  detail  microphysical  and  chemical  processes 
associated  with  PSCs.  This  model  describes  the  evolution  of 
the  microphysical,  thermodynamical,  and  chemical  properties 
of  an  air  parcel  following  an  isentropic  trajectory.  The  PSC 
code  and  the  photochemical  code  are  coupled  via  heterogene¬ 
ous  chemistry.  In  this  section,  the  different  packages  of  the 
model  are  described. 

2.1.1.  Trajectory  model.  MiPLaSMO  uses  isentropic 
trajectories  as  input.  Trajectories  are  calculated  [Knudsen  and 
Carver,  1994]  using  European  Centre  for  Medium-Range 
Weather  Forecast  (ECMWF)  meteorological  data  with  a 
spatial  resolution  of  1.125°  X  1.125°.  They  are  available  every 
6  hours.  Trajectory  data  are  calculated  every  2  hours. 
Mesoscale  perturbations  of  temperature  and  pressure  can  be 
prescribed  along  the  trajectory. 

2.1.2.  Microphysical  code.  The  microphysical  code  was 
initially  developed  at  the  Danish  Meteorological  Institute  by 
N.  Larsen  [Larsen  et  al.,  1997].  It  describes  the  time  evolution 
of  several  kinds  of  particles;  PSC  la  (supposed  to  be  NAT), 


Figure  1.  Scheme  of  PSC  formation  used  in  the  model  as  a 
function  of  the  temperature.  Typical  values  of  TNat,  Ticc  and 
7meit  are  indicated. 


PSC  II,  solid  sulfuric  aerosols  supposed  to  be  SAT,  and  liquid 
particles.  All  the  particles  are  distributed  in  50  size  bins  from 
0.001  pm  to  100  pm.  Figure  1  presents  the  pathways  used  in 
this  code  for  particle  evolution,  as  a  function  of  temperature. 

As  the  temperature  decreases,  the  liquid  binary  aerosols 
grow  as  STS  by  uptake  of  HNO3  and  H2O  from  the  gas  phase. 
PSC  II  particles  are  the  first  solid  particles  to  form  from 
freezing  of  liquid  droplets.  The  process  begins  at  temperatures 
less  than  Tlcc  -1  K.  According  to  Koop  et  al.  [1997],  particles 
can  remain  liquid  until  the  Tm  -3  K  threshold  is  reached. 
During  this  freezing  process,  the  condensed  HNO3  and  H2SO4 
are  conserved  through  the  formation  of  an  inner  NAT  shell 
surrounding  a  SAT  core  in  the  ice  particle.  All  the  sulfuric 
acid  is  in  the  condensed  phase.  As  the  temperature  warms  up, 
the  ice  particles  evaporate,  revealing  the  NAT  core  and  NAT 
can  evaporates  to  leave  SAT  particles.  Size  evolution  of 
liquid,  NAT,  and  ice  particles  due  to  the 
condensation/evaporation  mechanism  is  taken  into  account. 
Heterogeneous  nucleation  of  PSC  II  on  PSC  la  is  also 
considered  when  T-lce  is  reached.  In  the  model  we  assume  that 
SAT  can  return  to  liquid  particles  in  two  cases:  when 
temperatures  pass  beyond  the  melting  temperature  of  SAT,  or 
when  temperatures  go  below  the  temperature  of  SAT 
deliquescence  [Koop  and  Carslaw,  1996]  if  no  other  solid 
particles  exist. 

2.1.3.  Photochemical  code.  The  chemical  code  describes 
the  evolution  of  41  species  of  the  NOy,  HO*,  Cly,  O*,  and  Bry 
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families.  It  takes  into  account  112  homogeneous  reactions  and 
five  heterogeneous  reactions  on  clouds  or  aerosols.  Those 
heterogeneous  reactions  are 

(Rl)  HC1+C10N02->C12  +  HN03 

(R2)  H  20 + C10N02  ->  HOC1  +  HN03 

(R3)  N205+HC1->C1N02+HN03 

(R4)  N205  +  H20  -»  2  HN03 

(R5)  H0C1  +  HC1->C12+H20 

First-order  rate  constants  k  (s'1)  for  heterogeneous  reactions 
depend  on  type  and  surface  of  particles.  They  are  calculated  as 
follows: 

k  =  ySv  (1) 

where  S  is  the  total  area  of  a  particle  type  considered  for  the 
considered  particle  type,  v  the  mean  thermal  velocity,  and  y 
the  reaction  probability.  Reaction  probability  are  taken  from 
De  More  et  al.  [1997],  Carslaw  and  Peter  [1997],  and 
Ravishankara  and  Hanson  [1996]. 

The  numerical  scheme  used  for  the  chemical  code  is  semi- 
implicit  symmetric  (SIS).  The  nonlinear  equations  of  chemical 
species  evolution  are  converted  into  a  linear  system  of 
equations.  A  description  of  this  method  can  be  found  in  the 
work  by  Ramaroson  et  al.  [1992], 

2.2.  Mountain-Wave  Calculation 

Coolings  induced  by  mountain  waves  are  calculated  off¬ 
line  by  the  Naval  Research  Laboratory  Mountain  Wave 
Forecast  Model  (NRL/MWFM).  It  is  described  by  Bacmeister 
et  al.  [1994].  Extensions  of  the  model  to  calculate  temperature 
variability  are  described  by  Eckermann  and  Bacmeister 
[1998].  The  model  relies  on  a  database  of  dominant 
topographic  mountain  ranges  with  estimation  of  their  height, 
width,  and  orientation.  Vertical  motions  induced  by  the 
mountain  ridge  are  calculated  at  each  standard  pressure  level 
(including  50,  30,  and  10  hPa)  and  every  6  hours  by 
NRL/MWFM.  The  temperature  perturbation  is  then  deduced 
using  adiabatic  parcel  thermodynamics,  since  stratosphere 
mountain  wave  advection  is  adiabatic  to  a  good 
approximation  [Eckermann  et  al.,  1998].  The  temperature 
perturbation  introduced  in  MiPLaSMO  is  deduced  from 
cooling  pathways  previously  calculated  by  NRL/MWFM. 
Interpolations  are  made  to  the  pressures  and  times  of  interest. 
As  in  Tsias  et  al.  [1997],  the  mesoscale  temperature  variation 
is  assumed  sinusoidal.  The  amplitude  of  the  sinusoid  varies 
with  respect  to  the  maximum  cooling  computed  by  the 
mountain  wave  model.  The  period  of  the  oscillation  is 
estimated  from  the  mean  time  between  two  consecutive 
cooling  events  as  forecast  by  the  NRL/MWFM  model. 

3.  Instruments 

In  the  frame  of  the  validation  campaign  of  the  Improved 
Limb  Atmospheric  Spectrometer  (ILAS)  instrument  on  board 
the  ADEOS  satellite  [Sasano  et  al.,  1999],  chemical  species 
and  aerosol  measurements  were  performed  in  early  1997  with 
balloon-borne  instruments  from  Kiruna,  Sweden.  Here  we 
present  the  OPC,  AMON,  and  LPMA  instruments  which  have 
performed  measurements  on  February  25  and  26. 

3.1.  Optical  Particle  Counter 

The  OPC  used  for  in  situ  measurements  from  the  balloon- 
borne  gondola  measures  the  number  concentration  of 


condensation  nuclei  ( r  >  0.01  pm)  and  of  optically  detectable 
aerosol  (0.15  <  r  <  10.0  pm).  The  aerosol  counter,  following  an 
original  design  by  Rosen  [1964],  counts  and  sizes  aerosol 
particles  drawn  into  a  scattering  region.  The  size  is 
determined  from  measurements  of  the  intensity  of  scattered 
white  light  at  40°  from  the  forward  direction  using  Mie  theory 
and  assuming  spherical  particles  with  an  index  of  refraction  of 
1.45  [Hofmann  and  Deshler,  1991]. 

3.2.  AMON  Instrument 

AMON  is  a  balloon-borne  UV-visible  spectrometer 
devoted  to  nighttime  measurements  of  the  stratospheric 
species  O3,  NO2,  and  NO3  [Naudet  et  al.,  1994;  Renard  et  al., 
1996],  Aerosol  extinction  can  also  be  retrieved  as  a  function 
of  the  wavelength.  Observations  are  performed  using  the 
stellar  occultation  method,  which  consists  in  analyzing  the 
modifications  of  the  spectrum  of  a  setting  star  induced  by 
absorbing  atmospheric  species.  A  reference  spectrum  is 
obtained  when  the  star  is  above  the  balloon  horizon  and  the 
balloon  is  afloat.  The  transmission  spectra  are  obtained  by 
dividing  the  observed  spectra  by  the  reference  spectrum.  The 
retrieval  of  the  slant  column  densities  is  performed  using  a 
least  squares  fit  between  the  observed  transmission  spectra 
and  the  spectra  calculated  using  the  known  cross  sections  of 
the  absorbing  species.  The  residual  in  the  transmission 
spectra,  after  removing  the  contribution  of  Rayleigh  scattering 
and  the  contribution  of  molecular  absorption,  is  assigned  to 
aerosols.  The  vertical  profiles  of  the  species  and  of  the  aerosol 
extinction  are  obtained  after  inversion  performed  using  the 
onion  peeling  technique  [Renard  et  al.,  1996]. 

3.3.  LPMA  Instrument 

The  LPMA  is  a  remote  sensing  infrared  Fourier  transform 
instrument  operating  in  absorption  against  the  Sun  [Camy- 
Peyret,  1995].  With  its  high  spectral  resolution  and 
sensitivity,  the  retrieval  of  vertical  profiles  of  trace  species 
having  stratospheric  mixing  ratios  as  small  as  0.1  ppbv  is 
possible.  A  global  fit  algorithm  [Carlotti,  1988]  associated 
with  an  efficient  minimization  algorithm  of  the  Levenberg- 
Marquardt  type  [Press  et  al.,  1992]  is  used  for  the  retrieval. 
This  retrieval  algorithm  [Payan  et  al.,  1998]  allowed  us  to 
retrieve  vertical  profiles  of  HC1,  HNO3,  and  several  other 
stratospheric  species  from  occultation  data. 

4.  Study  of  February  25,  1997 

4.1.  OPC  Measurements 

Observations  were  performed  from  Kiruna  (67.9°N/22.1°E) 
on  February  25,  1997.  The  OPC  flight  started  at  0918  UT  and 
ended  at  1051  UT.  Data  were  recorded  from  the  ground  to  30 
km.  The  aerosol  concentration  profiles  for  different  size  classes 
are  shown  in  Figure  2.  A  layer  of  large  particles  was  observed 
between  20  and  22  km. 

4.2.  Initialization  of  MiPLaSMO 

For  this  case,  the  initial  distribution  of  particles,  typical  of 
background  aerosol,  is  deduced  from  OPC  measurements.  The 
distribution  at  the  desired  altitude  is  obtained  by 
interpolation  between  two  typical  distributions  of 
background  aerosols  observed  below  and  above  this  altitude. 
The  distribution  measured  on  February  23,  1997  (T.  Deshler, 
unpublished  ILAS  database,  1997)  has  been  used.  These 
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Figure  2.  Vertical  profiles  of  aerosol  concentration  for  particles  between  0.15  pm  and  10  pm  from  the  in  situ 
OPC  measurements  on  February  25,  1997. 


measurements  have  been  preferred  to  those  of  February  25 
because  the  disturbed  layer  involving  large  particles  is 
thinner.  This  allows  a  better  precision  in  the  interpolation  of 
aerosol  parameters.  At  initialization  of  the  model,  particles  are 
assumed  to  be  liquid.  Chemical  species  initialization  is  taken 
from  REPROBUS  3D  CTM  [Le/evre  et  al.,  1998]  simulations 
of  the  whole  winter  1996-1997. 

4.3.  Synoptic  Situation 

A  20-day  isentropic  back  trajectory  was  calculated  to 
obtain  the  temperature  history  at  the  500  K  level, 
corresponding  to  21.5  km  (35  hPa),  where  the  largest  particles 
are  observed  (see  Figure  3).  At  the  location  of  the 
measurement,  the  synoptic  temperature  (192  K)  is  5  K  over  the 
frost  point,  and  thus,  too  warm  for  a  PSC  formation.  Moreover, 
the  coldest  temperature  reached  along  this  trajectory  was  3  K 
above  the  ice  frost  point  (120  hours  before  the  measurement). 
As  a  consequence,  the  model  does  not  calculate  any  solid 
PSCs,  using  synoptic-scale  temperatures. 

4.4.  Mesoscale  Interpretation 

A  preliminary  study  of  this  case  was  made  by  Huret  et  al. 
[1998],  considering  a  theoretical  mesoscale  perturbation.  The 
authors  have  shown  that  the  hypothesis  of  a  mountain  wave 
event  could  be  an  explanation  to  interpret  the  measurements  of 
February  25,  1997.  Here  we  use  mountain-induced  wave  (lee 
wave)  perturbations  predicted  by  the  NRL/MWFM  model 
described  above  to  tentatively  explain  these  observations. 


At  500  K  the  OPC  measurements  were  expressed  in  terms  of 
concentration  as  a  function  of  the  radius.  The  back  trajectory 
for  February  25  crossed  the  Scandinavian  mountain  ranges 
prior  to  its  arrival  at  the  sampling  location  (not  shown).  Two 


Figure  3.  Synoptic  temperature  time  evolution  for  the  20-day 
back  trajectory  (heavy  solid  line)  at  the  500  K  level  for 
February  25,  1997.  Melting  temperature  of  SAT  (7mei,). 
temperature  of  NAT  condensation  (7nat),  temperature  of  SAT 
deliquescence  (Tdeiiq),  and  ice  frost  point  (Tice)  are  also  shown; 
t  =  0  corresponds  to  the  time  of  the  measurements. 
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Figure  4.  February  25  back  trajectory,  at  500  K.  (a)  Mesoscale 
temperature  perturbation  1  hour  before  the  OPC  measurements. 
Tmeit,  Tnat,  and  T\oc  are  presented  (see  Figure  3  caption). (b) 
Surface  area  time  evolution  of  the  particle  surface  area  during 
the  lee  wave  event  shown  in  Figure  4a.).  Liquid  particle, 
NAT,  and  PSC  II  total  surface  area  are  presented. 

hours  before  the  OPC  measurements,  NRL/MWFM 
calculations  show  that  the  air  parcel  temperature  had  been 
influenced  by  lee  waves,  generated  by  flow  over  the 
Scandinavian  mountains.  The  temperature  evolution  during 
the  last  hour  of  the  trajectory  is  shown  in  Figure  4a.  After  a 
succession  of  cooling-warming  cycles  of  5  K  amplitude,  the 
maximum  cooling  induced  by  the  wave  is  10  K,  which  occurs 
a  few  minutes  before  the  measurements,  and  T[Ce  -1  K  is 
reached.  Figure  4b  shows  the  evolution  of  particle  surface 
areas  during  the  lee  wave  event.  When  T-,Ce  -1  K  is  reached,  a 
small  amount  of  the  largest  STS  particles  freezes  into  PSC  II. 
Since  the  7jce  -3  K  threshold  is  not  reached,  some  particles 
remain  liquid.  Ice  particles  grow  by  condensation  of  water. 
Then,  a  warming  occurs  which  makes  the  ice  evaporate, 
leaving  NAT  particles.  A  cooling  occurs  again,  but  the 
temperature  stays  too  high  to  make  NAT  grow  by 
condensation  of  H2O  and  FINO3.  SAT  begin  to  appear  because 
of  NAT  evaporation  due  to  warm  temperatures.  Finally,  the 
distribution  computed  at  the  location  of  the  measurement  is 
composed  of  liquid  particles,  NAT,  and  low  concentration  of 
SAT  particles.  The  overall  size  distribution  resulting  from 
model  calculations  (Figure  5)  presents  a  bimodal  structure. 
This  bimodal  structure  is  in  good  agreement  with  the 
measured  spectrum  of  OPC.  The  first  mode  centered  on  a  radius 
of  0.013  pm  is  mainly  composed  of  liquid  particles  (with  a 


very  small  contribution  from  SAT  particles)  and  the  second 
mode  centered  on  1  p.m  is  due  to  NAT  particles. 

5.  Study  of  February  26,  1997 

5.1.  Measurements  by  AMON  and  LPMA 

Observations  were  also  performed  from  Kiruna,  on  February 
26,  1997.  The  AMON  observations  started  at  2105  UT  and 
ended  at  221 1  UT  with  a  balloon  float  altitude  of  31.5  km.  The 
star  used  for  the  occultation  was  Rigel,  and  100  spectra  were 
recorded.  The  latitude  and  longitude  were  in  the  67.7°-67.0°N 
and  27.3°-18.5°E  ranges,  respectively,  with  tangent  point 
altitudes  decreasing  from  31.4  to  11.3  km.  The  LPMA 
instrument  performed  a  flight  from  Esrange  (68°N/21°E),  just 
before  the  flight  of  AMON.  Infrared  spectra  were  recorded  at 
sunset  during  occultation  from  1500  UT  to  1540  UT.  The 
latitude  and  longitude  ranges  for  the  LPMA  occultation 
measurements  were  67.85°-66.87°N  and  26.77°-18.78°E, 
respectively,  and  so  were  very  close  to  those  of  AMON. 

5.1.1.  Aerosol  measurements  by  AMON.  Vertical  profiles 
of  aerosol  extinction  coefficient  obtained  every  25  nm  are 
shown  in  Plate  1.  Almost  all  the  wavelength  dependence  of 
the  aerosol  extinction  coefficients  at  different  altitudes  can  be 
fitted  using  a  Mie  theory  algorithm  including  a  log-log  size 
distribution.  The  comparison  of  the  aerosols  parameters 
derived  from  AMON  with  those  of  the  OPC  obtained  the  day 
before  shows  that  there  is  a  very  good  agreement  between  the 
two  measurements  in  terms  of  mean  radius.  In  the  altitude 
range  where  there  is  good  agreement,  the  retrieved  parameters 
correspond  to  background  aerosols. 

Significant  discrepancies  from  Mie  theory  occur  in  a  layer 
between  21  and  23  km  in  the  AMON  profile.  In  this  layer  the 
wavelength-dependent  extinction  is  maximum  for  the  larger 
wavelength  and  cannot  be  fitted  whatever  the  refraction  index 
and  size  distributions  used.  This  result  means  that  the 
hypothesis  of  spherical  particles  used  for  Mie  scattering  is  not 
valid.  Besides,  strong  chromatic  scintillation  was  observed 


radius  (nm) 

Figure  5.  February  25,  1997:  size  distribution  measured  by 
the  OPC  at  the  500  K  level  (stars)  and  computed  by 
MiPLaSMO  (solid  line).  Size  distribution  at  initialization 
used  for  the  simulation  is  also  presented  (dashed  line). 
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Plate  1.  February  26,  1997:  vertical  profiles  of  the  aerosol  extinction  coefficient  obtained  with  AMON  in  the 
near  UV  and  visible  domain  (from  375  nm  to  700  nm). 


by  AMON  in  the  layer  centered  around  22  km.  These  facts 
indicate  that  frozen  aerosols  are  probably  present  in  the 
corresponding  layer. 

5.1.2.  Chemical  trace  species  measurements  of  LPMA. 

At  the  level  of  the  layer  of  solid  particles  measured  by  AMON, 
the  LPMA  infrared  remote  sensing  measurements  of  HN03  and 
HC1  were  11.3±0.6  ppbv  and  0.24±0.11  ppbv,  respectively. 

5.2.  Synoptic  Situation 

The  model  was  initialized  in  the  same  way  as  described  in 
section  4.2.  We  calculate  the  temperature  history  of  the  air 
parcel  from  a  20-day  back  trajectory  at  the  520  K  level  (Figure 
6),  corresponding  to  21.6  km  (33  hPa),  where  the  extinction 
coefficient  for  the  larger  wavelength  is  maximum.  The 
temperature  at  the  location  of  the  measurement  is  197  K,  which 
is  1 1  K  above  7jce  and  more  than  3  K  above  7nat.  This  is  too 
warm  to  associate  this  solid  layer  with  the  presence  of  a  PSC. 
Furthermore,  the  coldest  temperature  encountered  along  the 
trajectory  is  more  than  4  K  above  the  ice  frost  point  (210 
hours  before  the  measurement),  so  we  conclude  from  the 
synoptic  temperature  history  that  no  PSC  should  have  formed 
because  of  the  warm  temperatures.  Furthermore,  no  solid 
particles  are  predicted  by  the  model  from  the  synoptic 
temperature  history. 

5.3.  Mesoscale  Interpretation 

Mountain-wave  model  calculations  were  performed  for  the 
two  times  when  the  sampled  air  parcel  crossed  mountain 
ranges.  On  February  26  over  the  Scandinavian  mountains,  the 
maximum  perturbation  amplitude  encountered  by  the  air  parcel 
is  5  K.  No  solid  particles  can  form,  since  the  lowest 
temperatures  are  still  5  K.  above  7jce.  However,  92  hours 
earlier  (February  23),  the  trajectories  indicate  that  this  air 
parcel  was  again  over  the  Scandinavian  mountains  and  the 
predicted  lee  wave  activity  was  much  stronger.  Figures  7a  and 
7b  show  the  mesoscale  temperatures  (computed  as  before)  and 
the  evolution  of  particle  surface  areas  during  this  period, 
respectively,  92  hours  prior  to  sampling  by  AMON.  A  13  K 
strong  cooling  was  first  predicted,  which  exceeds  the 


supercooling  threshold,  making  all  the  PSC  lb  particles  freeze 
into  ice.  Then  the  wanning  occurs,  whereupon  the  ice 
evaporates,  and  NAT  begins  to  appear.  NAT  also  evaporates 
in  this  wanning,  since  the  temperature  is  above  7nat-  The  air 
cools  again  to  188  K,  2  K  above  the  frost  point,  ruling  out  the 
possibility  of  ice  formation  from  heterogeneous  nucleation  on 
remaining  NAT.  The  remaining  NAT  particles  evaporate, 
leaving  SAT.  The  model  assumes  that  SAT  can  remain  frozen 
until  it  reaches  either  the  SAT  melting  temperature  or  SAT 
deliquescence  temperature.  Up  to  the  location  of  the 
measurements,  those  thresholds  are  never  reached  at  synoptic 
scale  (Figure  6).  Moreover,  when  the  air  parcel  crosses 
Scandinavia  just  before  the  AMON  measurements,  the 
aforementioned  mesoscale  processes  induce  cooling  and 
wanning  events,  but  the  maximum  cooling  induced  is  more 
than  2  K  above  the  SAT  deliquescence  temperature.  The 
maximum  warming  events  are  9  K  below  the  SAT  melting 


Figure  6.  Synoptic  temperature  time  evolution  for  the  20-day 
back  trajectory  (heavy  solid  line)  at  the  520  K  level  for 
February  26,  1997.  See  Figure  3  caption  for  the  temperature 
labels. 


RIVIERE  ET  AL.:  LEE  WAVES  AND  FORMATION  OF  SOLID  PSCs 


6851 


I . . . I 

-92  -91  -90  -89 

time  (h) 


Figure  7.  Same  as  in  Figure  4  but  for  February  26,  1997  at  520 
K,  92  hours  before  the  AMON  measurements:  (a)  temperature 
evolution  and  (b)  time  evolution  of  particle  surface  area. 

temperature,  so  SAT  particles  should  remain  solid.  When  the 
air  parcel  is  sampled  by  the  AMON  instrument,  SAT  is  the 
only  kind  of  particle  predicted  by  the  model. 

Regarding  the  chemical  composition  of  the  atmosphere,  the 
MiPLaSMO  microphysical  results  are  compatible  with  the 
high  HNO3  concentrations  measured  by  LPMA  [Payan  et  al., 
1999]  at  this  altitude,  ruling  out  the  presence  of  a  PSC  la.  The 
nitric  acid  computed  by  the  model  above  Kiruna  is  11.3  ppbv, 
while  the  amount  measured  by  LPMA  was  between  11.3±0.6 
ppbv  at  this  altitude  (Figure  8).  If  a  PSC  la  had  been  observed, 
the  amount  of  HNO3  would  have  been  lower,  due  to 
condensation  of  nitric  acid  and  water  from  the  gas  phase.  The 
HNO3  amount  computed  over  Kiruna  is  more  than  1  ppbv 
greater  than  the  amount  used  as  initialization.  This  is  due  to 
nitric  acid  production  by  heterogeneous  reactions  (Rl),  (R2), 
(R3),  and  (R4)  on  cloud  particles. 

The  model  interpretation  of  AMON  solid  particle 
measurements  is  then  related  to  a  previous  mesoscale  PSC 
event,  which  later  evaporation  led  to  the  release  of  frozen 
sulfuric  aerosol  particles. 

5.4.  Implications  for  Chlorine  Activation 

In  order  to  evaluate  the  impact  of  lee  wave  activity  on 
chlorine  activation,  model  runs  were  performed  with  and 
without  lee  waves  for  the  case  of  February  26,  1997.  Time 
evolution  of  HNO3,  HC1,  CIONO2,  and  CIO  are  shown  in 
Figure  8  for  cases  where  lee  waves  were  taken  into  account  or 
excluded.  The  initial  value  (REPROBUS  simulations)  of  CIO* 


(i.e.,  CIO  +  2  CI2O2)  is  equal  to  0.2  ppbv,  indicating  a  small 
chlorine  activation  of  the  air  parcel.  Around  210  hours  before 
the  measurements,  the  synoptic  temperature  was  3-4  K  above 
the  ice  frost  point  (Figure  6),  triggering  growth  of  liquid 
particles  by  absorption  of  nitric  acid  and  water.  This 
corresponds  to  an  increase  in  the  chlorine  activation  as  seen 
by  an  increase  of  the  calculated  CIO*  mixing  ratio  at  night  and 
a  calculated  decrease  of  HC1  of  0.36  ppbv.  This  period 
corresponds  to  a  slight  increase  of  nitric  acid  of  0.38  ppbv. 
Indeed,  the  cumulative  effects  of  large  particle  size  with  low 
temperature  and  the  liquid  phase  of  the  particles  induce  rather 
large  rate  constants  of  heterogeneous  reactions  and  thus, 
additional  chlorine  activation.  As  the  first  lee  wave  event 
occurs  over  the  Scandinavian  mountains  (at  around  /  =  -91 
hours),  the  air  parcel  is  already  activated  with  CIO  nighttime 
values  of  0.18  ppbv  and  CIO*  values  of  0.68  ppbv.  After  this 
event,  differences  of  more  than  a  factor  of  2  appear  in  daytime 
CIO  concentration  between  the  cases  with  and  without  lee 
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Figure  8.  Case  study  for  February  26,  1997,  at  520  K.  Time 
evolution  for  HNO3,  HC1,  CIONO2,  and  CIO  mixing  ratios, 
taking  into  account  mountain-wave  perturbations  starting  at 
about  t  =  -  91  hours  (dotted  gray  line),  and  excluding  them 
(solid  line).  HNO3  and  HC1  mixing  ratios  obtained  by 
REPROBUS  simulations  (triangle)  and  LPMA  measurements 
(gray  circle)  at  21  km  are  also  shown. 
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waves.  HC1  is  almost  completely  consumed  (0.1  ppbv 
remaining)  by  heterogereous  reactions,  and  CIONO2  decreases 
by  1.25  ppbv  for  the  same  reason.  The  CIONO2  loss  is  in 
correlation  with  the  nitric  acid  increase  of  the  same  amount 
due  to  reactions  (Rl)  and  (R2).  The  contribution  of  reactions 
(R3)  and  (R4)  for  nitric  acid  production  is  negligible  because 
of  the  very  small  amount  of  N2O5  in  our  conditions.  It  has  been 
verified  that  for  10-days  integration  period  using  synoptic 
temperature  history  that  REPROBUS  and  MiPLaSMO  results 
are  in  rather  good  agreement  (see  Figure  8  and  Table  1).  When 
lee  waves  are  introduced  into  the  model,  the  MiPLaSMO 
values  (at  21.6  km  at  the  time  of  the  AMON  measurements)  of 
HNO3  and  HC1  (11.3  ppbv  and  0.30  ppbv,  respectively)  are  in 
accordance  with  the  measurements  of  LPMA  on  February  26, 
1997,  at  21  km  (11.3±0.6  ppbv  and  0.24±0.11  ppbv, 
respectively).  Comparisons  between  MiPLaSMO  results  and 
LPMA  measurements  are  summarized  in  Table  2.  Differences 
between  LPMA  measurements  and  REPROBUS  simulations 
could  be  explained  by  the  fact  that  REPROBUS  is  not 
accounting  for  mesoscale  temperature  fluctuations. 

6.  Conclusion 

The  aim  of  the  present  study  was  to  interpret  two  sets  of 
measurements  performed  in  February  1997,  which  cannot  be 
explained  by  synoptic-scale  temperature  histories.  For  this 
purpose,  we  have  used  a  Lagrangian  chemical  and 
microphysical  model  (MiPLaSMO)  and  a  mountain  wave 
model  (NRL/MWFM).  We  have  shown  that,  introducing 
mesoscale  temperature  perturbations  from  NRL/MWFM, 
initialization  from  the  REPROBUS  model  (chemistry)  and 
OPC  measurements  (aerosols),  with  realistic  hypotheses 
about  the  supercooling  threshold  for  PSCs  freezing  (i.e.,  3  K 
below  the  ice  frost  point),  the  MiPLaSMO  Lagrangian  model 
could  produce  results  in  agreement  with  the  layer  of  large 
particles  observed  by  the  optical  particle  counter  on  February 
25,  1997.  This  layer  can  be  interpreted  as  composed  of  a  small 
number  of  NAT  particles  (PSC  la)  mixed  with  liquid  particles. 
The  maximum  temperature  amplitude  of  the  lee  waves  in  this 
simulation  was  10  K. 

On  February  26,  the  observation  of  solid  particles  can  be 
interpreted  with  MiPLaSMO  and  the  mountain-wave  model. 
Although  the  mountain-wave  model  is  considered  to 
overestimate  the  wave  amplitude  [Carslaw  et  al.,  1999],  the 
13  K  lee  wave  temperature  amplitude  used  in  our  model  is 
nevertheless  still  compatible  with  the  observational  case 
studied  by  Carslaw  et  al.  [1998a].  What  is  observed  could  be 
solid  sulfuric  acid  aerosol,  evidence  of  a  previous  PSC  event. 


Table  1.  Comparison  Between  REPROBUS  at  Two  Standard 
Pressure  Levels  (21  km  and  21.9  km)  and  MiPLaSMO  (21.6 
km)  Model  Results  for  HNO3  and  HC1. 


HNO3  (ppbv) 

HC1  (ppbv) 

REPROBUS  (LPMA  location) 

(21  km) 

MiPLaSMO  (21  6  km) 

9.7 

1.00 

without  lee  waves 
(AMON  location) 

10  3 

0  96 

REPROBUS  (LPMA  location) 
(21.9  km) 

9.6 

1.30 

Temperature  and  pressure  perturbations  are  not  taken  into  account 


Table  2.  Comparison  Between  LPMA  (1530  UT)  or  AMON 
(2200  UT)  Measurements  at  21  and  22  km,  and  MiPLaSMO 
Model  Results  at  21.6  km  for  HNO3  and  HC1. 


HNO3  (ppbv) 

HC1  (ppbv) 

LPMA  (21  km) 

1 1  3  ±0.6 

0.24±0  1 1 

MiPLaSMO  (21.6  km) 

with  lee  waves 

11.3 

0.30 

(AMON  location) 

LPMA  (22  km) 

11.0±0.6 

0.45±0.20 

Here,  mesoscale  perturbations  due  to  lee  wave  are  introduced  in 


MiPLaSMO. 


This  PSC  would  have  been  generated  by  wave  perturbations 
above  the  Scandinavian  mountains,  several  days  before  the 
measurements.  Solid  aerosol  calculations  are  qualitatively  in 
accordance  with  the  high  HNO3  concentration  and  the  low 
values  of  HC1  measured  on  February  26.  For  this  case,  we  have 
also  shown  that  even  if  the  air  parcel  was  already  chlorine 
activated  before  the  mesoscale  processes  occurred,  the  lee 
wave  activity  has  significantly  amplified  the  activation 
process  and  CIO  concentration  is  more  than  2  times  greater 
than  without  lee  waves  in  daytime  conditions.  This  result  is 
in  agreement  with  the  conclusion  of  Carslaw  et  al.  [1998b], 
Chemical  results  of  MiPLaSMO  are  in  accordance  with  HNO3 
and  HC1  LPMA  measurements.  Furthermore,  mesoscale 
processes  can  be  an  explanation  for  the  differences  observed 
by  Payan  et  al.  [1999]  between  the  REPROBUS  model  and 
LPMA  measurements. 
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